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Interference at a beam splitter reveals both 
classical and quantum properties of electromag- 
netic radiation. When two indistinguishable sin- 
gle photons impinge at the two inputs of a beam 
splitter they coalesce into a pair of photons ap- 
pearing in either one of its two outputs. This 
effect is due to the bosonic nature of photons and 
was first experimentally observed by Hong, Ou, 
and Mandel (HOM) [1]. Here, we present the 
observation of the HOM effect with two indepen- 
dent single-photon sources in the microwave fre- 
quency domain. We probe the indistinguishabil- 
ity of single photons, created with a controllable 
delay, in time-resolved second-order cross- and 
auto-correlation function measurements. Using 
quadrature amplitude detection we are able to re- 
solve different photon numbers and detect coher- 
ence in and between the output arms. This mea- 
surement scheme allows us to observe the HOM 
effect and, in addition, to fully characterize the 
two-mode entanglement of the spatially separated 
beam splitter output modes. Our experiments 
constitute a first step towcu-ds using two-photon 
interference at microwave frequencies for quan- 
tum communication and information processing, 
e.g. for distributing entanglement between nodes 
of a quantum network [2, 3] and for linear optics 
quantum computation [4, 5]. 

Presently, HOM two-photon interference has been 
demonstrated exclusively using photons at optical or tele- 
com wavelengths. Experiments were performed with 
photons emitted from a single source using parametric 
downconversion [1], trapped ions [3], atoms [6], quantum 
dots [7] and single molecules [8]. The HOM effect has 
also been observed with two independent sources [9-15] 
realizing indistinguishable single-photon states which are 
required as a resource in quantum networks or linear op- 
tics quantum computation. Such experiments have also 
been performed using donor impurities as sources [16] 
including NV-centers in diamond [17, 18]. Furthermore, 
the HOM effect has been employed to create entangle- 
ment between ions [19] in spatially separated traps, and 
to realize a controUed-NOT gate in a small-scale photonic 
network [20]. 

Here, we demonstrate the HOM interference of two 



indistinguishable microwave photons emitted from in- 
dependent triggered sources realized in superconducting 
circuits. In each source A (B) a transmon-type qubit is 
strongly coupled with rate g/in = 169 (177) MHz to a 
transmission line resonator, sec Fig. la. To create a sin- 
gle photon using one of the sources, we coherently excite 
the qubit into a state a\g) + 13 \e) using a resonant mi- 
crowave piilsc. Tlicn, we swap the qubit state into the 
resonator mode A (B) by tuning the qubit transition fre- 
quency for half a vacuum-Rabi period into resonance with 
the resonator. This creates the state a ]0) + /? Jl) in the 
resonator, which for /3 = 1 corresponds to a single-photon 
Fock state [22]. The photon then decays exponentially 
with Lorentzian spectrum through the strongly coupled 
output port of the resonator into the input mode a' (6') 
of the beam splitter at a rate k/2tt = 4.1 (4.6) MHz, see 
Fig. 1. The two photons then interfere at the beam split- 
ter and are emitted into the output modes a and b, see 
Fig. la. Using the dispersive interaction between qubit 
and resonator, we tune the emission frequencies of the 
two sources to an identical value of = 7.2506 GHz. 
This is achieved by adjusting the qubit transition fre- 
quencies Va to 8.575 (8.970) GHz with a static magnetic 
flux applied to the SQUID-loop of each qubit. For our 
experiments, we sequentially create 20 single photons in 
each source at a rate 1/tr = 1/512 ns ~ 1.95 MHz in a 
sequence repeated every 12.5 fis. 

To probe the photon statistics in the beam splitter 
outpiit modes a, and h we use two spatially separated 
heterodyne detection channels [21, 23], see dashed box 
in Fig. lb. Each channel consists of a set of semicon- 
ductor linear amplifiers, a microwave frequency mixer 
for downconversion to 25 MHz, followed by an analog- 
to-digital converter. After further digital downconver- 
sion to (i.e. and filtering, we extract the two quadra- 
ture amplitudes Xa,/b{t) and Pa/b{t) corresponding to 
the complex envelope of the two time-dependent signals 
Sa/b{t) — Xa/b{t) + iPa/b{t) [24]. In contrast to many 
other HOM experiments in which photons in the beam 
splitter outputs are detected by single-photon counters, 
our measurement of the complex envelope is intrinsically 
photon-number resolving for averaged measurements and 
allows us to measure coherences of the electromagnetic 
field. The complex envelopes Sa/b{t) a-re used to com- 
pute the statistics required to extract all relevant quan- 
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Figure 1. Sample and schematic, a, False color micro- 
graph of the HOM device. A transmon qubit (blue, see inset) 
controlled by a flux (purple) and charge gate line (orange) 
through which current or microwave pulses are applied, re- 
spectively, is coupled to each coplanar waveguide resonator 
(green) with fundamental mode AiB). The coupling rate of 
the input port of the asymmetric resonator is by a factor of 
10^ smaller than at the output port. Each resonator output is 
coupled to the input mode a! (6') (inset) of a microwave beam 
splitter (dark red) with two output modes a (b). b. Schematic 
of the HOM experiment with two individual cavity QED sys- 
tems for photon generation, a beam splitter (BS), and linear 
amplification and heterodyne detection for correlation func- 
tion measurements and full quantum state analysis [21]. 



turn correlations measured here [21, 24]. To analyze the 
digital data, signal processing and statistical analysis are 
performed in real-time using field programmable gate ar- 
ray (FPGA) based electronics [22]. The noise added by 
the detection chain is fully characterized by measuring its 
statistical properties when the output modes a and b of 
the beam splitter are left in the vacuum state with both 
sources idle [21, 24]. To verify the single-photon charac- 
ter of each source individually, we have measured their 
second-order cross-correlation functions G)^^ (r). These 
display clear antibunching [22], see Fig. 2ab. 

To investigate two-photon quantum interference, we si- 
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Figure 2. Second-order correlation function measure- 
ments. Cross-correlation G^^^(r) of individual single photon 
sources displaying single photon antibunching, source A in 
panel a and B in panel b. a-e, G^^^j (r) when operating both 
sources simultaneously with delay times St = 0,50,150 ns 
displaying the HOM effect and its dependence on the tem- 
poral photon overlap at the beam splitter, f, Second-order 
auto-correlation Gio* (t) of mode a for 5t — displaying two- 
photon coalescence. See main text for details and a discussion 
of the offsets indicated in the upper right corner of each panel. 
Blue dots are data and red lines are theory. 



multaneously generate two indistinguishable photons ide- 
ally realizing a two- mode entangled state (|20) + \02)) / \/2 
at the beam splitter outputs. The measured cross- 
correlation of the beam splitter output powers is observed 

(2) 

to vanish G\^{t) w for all r between —tr/2 and tr/2, 
see colored region in Fig. 2c. Therefore, we conclude that 
both microwave frequency photons coalesce at the beam 
splitter. This is the HOM effect with microwave photons. 
In our experiments the spatio-temporal coherence of the 
single-photon states is governed by resonator decay alone, 
and shows no significant additional dephasing. This is 
in stark contrast to many other experiments in which 
decoherence resulting in random frequency differences 
between interfering photons causes finite correlations at 
T-0 [6, 11, 13, 18]. At T^ntr (n = ±1, ±2, . . . , ±10), 
the peak at integer non-zero multiples of the photon gen- 
eration period reflects the product of the power in each 
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output (a^o) {h^b) which wc have normahzcd to one. 

All measured second-order correlation functions are 
normalized by a common scaling factor and a relative 
gain between the two amplification channels. A remain- 
ing small offset subtracted from the reconstructed corre- 
lation functions is indicated together with its standard 
deviation (extracted from multiple measurements) in the 
upper right corner of each panel in Fig. 2. The measured 
correlation functions are in good agreement with analyt- 
ical calculations [25] (solid lines in Fig. 2), taking into 
account the cavity decay rates k extracted from indepen- 
dent measurements and a fixed detection bandwidth of 
20 MHz chosen to reject experimental noise outside the 
desired band. 

To explore the level of indistinguishability between 

the two interfering photons we introduce a time delay 
6t on the order of the photon decay time 1/k ~ 37 ns. 

(2) 

For St = 50 ns the correlation function G)^^ (t) remains 
close to zero at r = 0, indicating the coalescence of 
those photons detected with vanishingly small time dif- 
ference [6, 25, 26], see Fig. 2d. Typically this effect is diffi- 
cult to observe with dcitcctors of insufficient bandwidth or 
sources with significant dephasing rates [13]. The small 
positive correlations observed at r ~ are due to the 
decreased temporal overlap of the single-photon spatio- 
temporal mode functions at the beam splitter, i.e. the in- 
creased distinguishability of the two photons. At r ~ n 
we observe broadened, lower amplitude correlations. All 
features are in agreement with theory (solid lines) [25]. 

For St = 150 ns (and St = 100 ns, not shown) the en- 
velopes of the two single-photon mode functions barely 
overlap at the beam splitter, resulting in fully distin- 
guishable single photons. At r = ±St and t = ntr±ST, 
see Fig. 2e, we observe positive correlations with an am- 
plitude 1 /4, which originate from single photons imping- 
ing on the beam splitter at different times and at different 
input arms, compare Fig. 2ab. At t = ntr the correla- 
tions between photons in the same beam splitter input 
arm sum up to 1/2, as expected. 

To clearly distinguish between single-photon anti- 
bunching and two-photon coalescence in time-resolved 
correlation function measurements, we have also mea- 
sured the second-order auto-correlation function G^a (t) 
of mode a. When operating only one single-photon source 

(2) 

Gaa (0) is expected to vanish. However, in the HOM con- 
figuration with St = we find (jia'*(0) = G^aa{ntr) = 
G\^{ntr) = 1, see Fig. 2f, as there is a 50% probability 
of detecting two photons in mode a. All measurements 
of Gaa{T) are in good agreement with calculations [25], 
both for St = Q (solid lines) and for St = 100 ns (not 
shown). Here, it is interesting to note that a second-order 
auto-correlation function is rarely directly measured, be- 
cause of the the lack of sufficiently fast single-photon de- 
tectors [27]. In contrast, using our heterodyne detection 
scheme we are capable of measuring Gia'(r) for multi- 



photon states. 

To distinguish between an equal mixture of the states 
]20) and |02), compatible with the observed correlations, 
and their coherent superposition ([20) + )02))/\/2, we 
fully characterize the two-photon states created in our 
HOM experiment by quantum state tomography. This 
allows us to probe the entanglement generated between 
the coalescing two-photon states in the two output ports 
of the beam-splitter. The created states are also referred 
to as NOON states. In contrast to the NOON states 
of propagating photons investigated here, NOON states 
have also been investigated in superconducting circuits 
with photons localized in resonators [28, 29]. 

To perform full quantum state tomography on prop- 
agating photons, we record four-dimensional histograms 
of the measured field quadratures Xa, Pa, Xi,, and Ph- 
This is an extension of the scheme discussed in [21] to 
two spatially separated modes. Prom these measure- 
ments we extract all moments of the two-mode field, 
i.e. expectation values of the form {{d^)"d"^{b^)''b'') with 
n, m,k,l G {0, 1, 2}. We observe that the first-order mo- 
ments and all other odd-order moments are zero since 
all single-photon Fock states are characterized by a fixed 
photon number and consequently a fully random phase, 
see Fig. 3a. Since each mode carries exactly one pho- 
ton on average, {d^d) and (b^b) are close to unity, while 
all other second-order moments vanish. The fourth-order 
moment {d^d¥b) is observed to be zero, while (a^a^ dd) 
and {b'b'^ bb) are unity, consistent with the coalescence of 
the two photons into either output. The above observa- 
tions are consistent with the ones based on the correla- 
tion function measurements in Fig. 2. Most importantly, 
the two-mode entanglement is indicated by the moment 
(ddww) which is close to 1, as expected. All measured 
moments of the two-mode entangled state created in our 
HOM experiment are in good agreement with the pre- 
dicted ones, sec wire frames in Fig. 3a. Note that mo- 
ments of order five and higher are all close to zero within 
their statistical errors. 

In addition, we have determined the most likely density 
matrix p characterizing the created two-mode entangled 
propagating photon state from the measured moments 
and their respective standard deviation following [21]. 
We have restricted the evaluation of moments to less than 
three photons per output mode, since we create no more 
than two single photons with our sources. Note that a 
related analysis has recently been performed in Ref. [30] . 
The real part of p is shown in Fig. 3b, all elements of the 
imaginary part of p are smaller than 0.02 (not shown). 
We extract a fidelity of the NOON type state of 84 % and 
a negativity of 0.39. 

Finally, to explore the interplay between single- and 
two-photon interference, we have performed experiments 
with modes A and B prepared in superpositions of 
and 1 photon Fock states, i.e. (|0) — i|l))/\/2 and 
(|0) -|- e"^ )l))/-\/2 with variable phase (p, ideally creat- 
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Figure 3. Full quantum state tomography of propagating photon states created in HOM experiments, a, Mea- 
sured moments {(a^)"a'"(b^)*6') with n,m,k,l G {0, 1, 2} up to fourth order and b, density matrix (real part) for (|20) + i02))/\/2 
with fidelity 84%. Colored bars are extracted from data, error bars indicate the standard deviation and wire frames represent 
ideal state, c, Measured moments for (V2|00) + 2 |01) + |20) + |02))/\/8 with corresponding state fidelity 84%. d, Power 
detected in the two output modes ((a^a), blue) and {{o'b), red) for operating only source B (line) or both sources (dots) 
creating superposition states with relative phase angle ip, see text for details. 



ing the state (^2 |00) - i{l - e''^) |10) + (1 + e"^) |01) + 
e"^(|20) + \02)))/V8 at the beam splitter output. With 
these input states, we have first measured the power in 
the beam splitter output modes a (blue line) and b (red 
line) when only operating source B and keeping source A 
idle, see Fig. 3d. In this case, we observe a power level 
corresponding to 1/4 of a single photon independent of 
the phase angle (p, as expected for an equal superposi- 
tion of |0) and |1) impinging on a balanced beam splitter. 
Operating both sources A and B, we observe a sinusoidal 
interference with phase ip of the two superposition states 
in the beam splitter output power of mode d (blue dots) 
and b (red dots), respectively. The sinusoidal oscillation 
is a result of the interference between one photon in ei- 
ther output port while the offset in power of 1/4 is the 
result of two-photon coalescence. 

For the phase angle « we have also performed 
full quantum state tomography, see Fig. 3c. As for the 
two-photon NOON state we observe antibunching, co- 
alescence, and entanglement in the moments (a^ab^b), 
{a^a^dd) and (b^b^bb), and (ddb^b^), respectively with 
close to expected amplitudes (wire frames). The inter- 
ference of the superposition states is revealed not only 
in the power {d^d) and (6^6) but also in the coherences 
with an unbalanced number of creation and annihilation 
operators. The corresponding state has a fidelity of 84 % 
with respect to the ideal one. 

Our results suggest that multiple on-demand single- 
photon sources emitting indistinguishable single photons 
could be used for creating nonlocal entanglement in quan- 



tum repeater or quantum communication applications 
based on microwave photons. 
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